o 
o 

(N 



> 
On 
l> 
O 

o 
t> 
o 

Or 



X 



CONTINUOUS SYMMETRIES OF THE LATTICE POTENTIAL 

KDV EQUATION 

DECIO LEVI AND MATTEO PETRERA 



Abstract. In this paper we present a set of results on the integration and on 
the symmetries of the lattice potential Korteweg-de Vries (lpKdV) equation. 
Using its associated spectral problem we construct the soliton solutions and the 
Lax technique enables us to provide infinite sequences of generalized symme- 
tries. Finally, using a discrete symmetry of the lpKdV equation, we construct 
a large class of non-autonomous generalized symmetries. 



1. Introduction 
The lattice version of the potential Korteweg-de Vries (lpKdV) equation 

(1) w t =w xxx + 3wl, 

is given by the nonlinear partial difference equation |17j : 

(2) ES> = (p - q + U n ,m+1 - U n +l,m) (p + £? - M„+l,m+l + U n , m ) - (p 2 - q 2 ) = 0. 
The above equation is probably the best known completely discrete equation which 
involves just four points which lay on two orthogonal infinite lattices and are the 
vertices of an elementary square - a quad-graph - (see Fig.l). Eq. ((2]) is one of the 
lattice equations on quad-graphs classified in [I], where the 31? consistency is used 
as a tool to establish its integrability. 



*n,m-(-l 



P 



Un+1,. 



Figure 1. An elementary square 



The lpKdV equation has been introduced for the first time by Hirota [8] in 1977 
and it is nothing else but the nonlinear superposition formula for the Korteweg-de 
Vries equation under disguise. A review of results about the lpKdV equation can be 
found in [T7] . 
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Integrable equations possess an infinite set of symmetries. Few of them arc point 
symmetries, i.e. symmetries whose infinitesimal generators depend just on the in- 
dependent and dependent variables, while an infinite denumerable number of them 
are generalized symmetries. The latter ones depend on the derivatives of the de- 
pendent variable with respect to the continuous independent variables and on a few 
lattice points if the independent variables are discrete. The presence of this infinite 
Lie algebra of symmetries and of the associated conserved quantities is one of the 
most important features of the integrability of a given nonlinear equation and it has 
been used with profit in the past to provide integrability tests for several partial 
differential equations in R 2 and differential-difference equations [TSJ [24] . 

To be able to introduce an integrability test based on symmetry one needs to 
understand the structure of the infinite dimensional symmetry algebra of integrable 
equations. In the case of completely discrete equations the situation is not as clear as 
for the differential-difference or the partial-differential case. A result in this direction 
has been obtained for the discrete-time Toda lattice [TT]. However, the Toda lattice 
is just an example and more examples are needed to get a sufficiently general idea 
of the possible structures which may appear. 

On the other hand, in [9], [10l [12] , the multiscale expansion technique [22] has been 
extended to the case of equations defined on a lattice. In [10] one of the authors has 
performed, as an example, the multiscale expansion of eq. ^) deriving a completely 
local discrete nonlinear Schrodingcr equation (dNLS). The integrability of this dNLS 
equation has been questioned by many people. So the construction of the symmetries 
of the lpKdV equation and their multiscale expansion should be a concrete tool to 
check its integrability. 

The present paper is devoted to the study of the lpKdV equation exactly with 
these aims. In Section [2] we review some known results on this lattice equation, 
while in Section [3] we present its inverse scattering transform, which turns out to 
be a slight generalization of the results presented by Boiti and collaborators on 
an asymmetric discretization of the continuous Schrodinger spectral problem [3]. 
Section [4] is devoted to the construction of Lie point and generalized symmetries. 
Finally, Section [5] contains some concluding remarks and some open problems. 

2. The lattice potential KdV equation 

Here we present some known results on the lpKdV equation ([2]) and its integra- 
bility. 

In eq. ([2]) u„ jm is the dynamical field variable, which we assume to be real and 
asymptotically bounded by a constant, defined at the site (to, n) £ I? while p, q are 
two nonzero real parameters which are related to the lattice steps and will go to 
infinity in the continuous limit so as to obtain from eq. ([2]) the continuous potential 
KdV equation U). 

Since we have two discrete independent variables, n and to, we perform the con- 
tinuous limit in two steps. Each step is achieved by shrinking the corresponding 
lattice step to zero and sending to infinity the number of lattice points. 

In the first step we define u n _ rn = Vk(r), where k = n + m and r = <5to, being 
S = p — q the lattice step in the m-direction. In the limit to — > oo, 5 — > 0, we obtain 
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the differential-difference equation 

(3) = 2p - 1. 

dr 2p-v k+1 +v k _ 1 

Defining 

q k = 2p- v k+2 +V k , 
we can rewrite eq. ^ as the differential-difference equation 

/ % dq k ( f 1 

(4) # = 2 / 



dr \qk-i qk+i, 

Eq. (H} has been obtained in ref. [5] as the simplest local evolution equation asso- 
ciated with the asymmetric discrete Schrodinger spectral problem 

(5) ip k+ 2 = qk tpk+i + A^fe, 

where A € C is a spectral parameter. 

Introducing the new field s k = (2p)/q k , eq. (jl]) reads 

(6) — = s fe (s fe+ i - Sfe-x), 

the so called discrete KdV equation [16 . Then the Miura transformation a k = 
s fe s fc-i: maps eq. ([5]) into the discrete Volterra equation 

(7) = a k {a k +i - ajfe-i), 

associated with the Schrodinger spectral problem 

ipk-i +a k ip k+1 = fitpk, 
where /i € C is a spectral parameter. 

The second continuous limit of eq. @ is performed by defining v k (r) = w(x,t), 
with x = 2 (fc + r/p) /p and t = 2 (fc/3 + r/p) /p 3 . If we carry out the limit p — > 
oo, fc — > oo, t — ► oo, in such a way that x and i remain finite, then eq. ([3]) is 
transformed into the potential KdV equation 

The integrability of the lpKdV equation j2]) is proven in ref. [17] by giving its 
Lax pair, an overdetermined system of matrix equations for the vector m = 
(V'n mQ 1 ), ipn m (h)) T , where ft, G C is a spectral parameter: 
(Rn} \& h — T h \b h 

\ u 1 n+l,m n,m n,m> 

n,m+l n,m n,m> 



(8b) , , = M' 1 *' 1 



with 



and 



P ~~ u n+l,m 1 

h 2 - p 2 + (p + u,^ m ) (p - u n+1 . m ) p + u ri 



M h = 

n.m 



q — Un,m+1 1 
h 2 - q 2 + (q + U n ,m) (Q - %,ro+l) <7 + " 
The consistency of eqs. ([8]) implies the discrete Lax equation 



n :in 



T h M h = M h T h 

n,m+l n.m n+l,m n,ml 



which plays the same role as the zero-curvature equation L it = [M , L ] in the con- 
tinuous case. The lpKdV equation © corresponds to an isospectral deformation of 
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the Lax pair ([8]), i.e. whenever h is an m and n-independent complex constant. We 
can rewrite eqs. ([8]) in scalar form in terms of ipn,m = V'i mW : 

(9a) ^n+1,m = (2p- U n+2i m + U n<m ) Ipn+l.m + (h 2 - p 2 ) 1p n ,m, 

(9b) 1pn, m +l = 1pn+l,m + (q~P + U n + i. m - M„, m+ l) i> n ,m- 

Taking into account eqs. (|9bf f2")l we can rewrite the m-evolution also as 

(10) 1pn, m +2 = {2q- U„, m+2 + U n<m ) 1p n ,m+l + (h 2 - q 2 ) V>n,m- 

Eq. (|10[) coincides with eq. (|9ap swapping the indices n <-> m and p <-» q. It is easy 
to see that also the lpKdV equation has the same discrete symmetry. 
Defining 

(11) q n = q n ,m =2p — U n +2,m + Un,m, 

choosing A = h 2 —p 2 and dropping the parametric dependence on m we see that eq. 
(f9a|) is equivalent to the discrete Schrodinger spectral problem ([5]). The m-evolution, 
provided by eq. (|9b[) . cannot be written in a simple way in terms of q n ^ m - 



3. The discrete spectral problem associated with the lpKdV equation 

Let us now study the direct and inverse problems associated with eq. ([5]). Our 
results are a generalization of those contained in [3|: in fact they reduce to them 
when p = 1. Henceforth we remand to [2] for most of the technical details. 

As the independent discrete variable m - the discrete time variable of the lpKdV 
equation - enters parametrically in the results of this section, it will not be explicitly 
written. We shall write it just when it is necessary, namely when we will study the 
m-evolution of the spectral data associated with the lpKdV equation. 

3.1. Direct problem. The solutions u nyTn of eq. ^ must go asymptotically to 
an arbitrary constant to be consistent with the difference equation ((2|). Then q n , 
defined in eq. (fTTj) , goes asymptotically to 2 p. 
Following [3] we define 

(12a) h = in, \ = -n 2 -p 2 , 

(12b) q n = r) n + 2p, ip n = (p + i«) n x„, 

where r\ n = —u n+ 2.m + i„,m is a new asymptotically vanishing potential. In terms 
of the new variables (fl2| the spectral problem ([5]) reads 

(13) (p + lK)Xn+2 - 2pXn+l + (p-i«)Xn = VnXn+1- 

The Jost solutions /i„ of the spectral problem fTB")) are defined in terms of the 
potential rj n and of the discrete complex exponential \{p + i/c)/(p — \n)] n through 
the following "discrete integral equations" 

i +oo / . \ j—n 

( 14a ) ti = 1 -2ik £ 1 



j=n+l 

(14b) ^=l + i~ K t 

3=— oo 



Vj-i Mj 
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The Jost solution fj,„ is an analytic function of n for lm(n) > and fi n for Im(ft) < 
with the boundary conditions 

(15a) lim fj+ = 1, Im(/c) > 0, 

n — *-}-qo 

(15b) lim = 1, Im(/c) < 0. 

n — > — oo 

For Im(«;) = we define the spectral data: 



3=-oo 
+°° / m i ; ~ \ J 



J=-oo 

Due to the analiticity property of the Jost solutions one can prove that a + (k) can 
be analytically extended to Im(ft) > and a~(n) to Im(«;) < 0. The inverse of 
the functions a ± (K), namely T ± (k) = [a play the role of the transmission 

coefficients and the functions R (k) = ^ (re) /c^ (/t) are the reflection coefficients. 
The poles of T ± (k) are related to the soliton solutions of the evolution equations 
associated with the spectral problem lfT5|) . Assuming that eqs. (|14l) are solvable and 
their solution is unique, we obtain for Im(re) = the following relations between the 
Jost solutions and the spectral data, 

Mn(«0 = ± ( K )^,T( K ) + ( , K ) b ± (K,)f^(-K). 



± 



y p + 1 K, 

Taking into account the asymptotic limits (| 1 5[) of the Jost functions, for Im(re) = 
we get 



(17) 



1 p + lK ' 



=poo. 



Let us assume that the function o + (k) has N simple zeros at n = kT, 1 < j < AT. In 
correspondence with these values of n the Jost function is asymptotically bounded 
and the N residues of the transmission coefficient T + (k), say {c'j'}jL i , are associated 
with its normalization coefficient. 

Henceforth we can associate with a given potential r\ n in a unique way the spectral 
data S[r) n ] = {R(k), T(k), {cj}j =1 }, obtained as a solution of the spectral problem 
(fT5|l . As the spectral problem (|13[) is a linear ordinary difference equation in n, the 
solution Xn is defined up to an arbitrary normalization function fl m (K). 

3.2. Inverse problem. Taking into account eqs. (fT5HT7j) and using the Cauchy- 
Green formula we can write the following integral equations for the two Jost solutions 



(i 
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(18a) 



N 

3 = 1 
1 



(p 



^n(- s )R( s ) (P 



1 s 



i \p + i s 



ds, Im(ft) < 0, 



N 



3 + 
K — Kj 
3 = 1 ■? 



V + 1K 



(18b) 
where 

(19) 



1 

2Tk 



iO 



p + is 



ds, Im(«;) > 0, 



E 



.+ 



2iK, , « — ' \ p - IK- 

J k— — oo \ 3 



vk-ivt(Kj)- 



Therefore, given the spectral data S [i] n ], we are able to reconstruct the Jost solutions 
and, in terms of them, the potential 

2 



N 



Vn 



(20) 



Mn+2 (-4) 



Mn+2,m( S ) 



p — i/d 



P + 1 S 
p — i S 



Mn,m( K j ) 



P + IK, 



P - 1 K 



Mr, 



p + 1 s 

p — is 



R+(s)ds. 



We remand to [3] for further details and for a study of the convergence of the sum 
appearing in the definition (fin)) of the coefficients Cj . 



3.3. Time evolution of spectral data. Taking into account the definitions (|T2|) , 
the m-evolution of Xn,m is obtained from eq. (|9b[) and is given by 

(21) Xn,m+1 = (p + i K) Xn+l,m + ~ P + "n+l.m ~ U n<m+1 ) Xn,m- 

As the lpKdV equation ^ is an isospectral deformation of eq. (fT3")) . from (|12ap « 
is an m-independent parameter. 

Defining Xn,m = ^mWrfm, where fi m (K) is a normalization function, and tak- 
ing into account the asymptotic behaviours (|15H17[) of the Jost functions, we find 
from eq. (l2Tj) the following time evolution of S[r/„ jm ] 



(22) T m+l ( K )=T m (K), 
Integrating eqs. (|22|) we get 

T m (re) = T(k), 



flm+l(t) 



R m (n). 



R(k), 



Rm(K) 

q + lK / 

namely the transmission coefficient is invariant under the evolution of the lpKdV 
equation while the reflection coefficient adquires an m-dependent discrete expo- 
nential factor. 
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3.4. Soliton solutions. The soliton solutions of eq. ([2|) are obtained by considering 
the reflectionless solution of the spectral problem (l9al [9bl IT2|) . They are obtained 
from eq. (|20|) setting = 0. The Jost solutions /4~ m (— i«t) are obtained 

solving eq. (|18b|) for k = — For AT = 1 they are given by 

_i_ \ iin — l 

p — IK 



while for N = 2 they read 



1 



2k+ Vp + i/^ 



where 



,(-4) 



Co 



1 


n 


c 2 


m 




i«2 


X>1,2 


4 - 


4 


\p + 




1 


H 




m 


fp- 


i 4 


©1,2 


4- 


K 2 




i 4 



P — 1 « 2 



Ci, 



2 \p + i4 



v - 



2k{ 

p — IK 



p + i K{ 
P 



44 4(4 + 4)2 

Taking into account the relation between rj n 
we get that the one soliton solution for the lpKdV equation 

q ( p— i K + \ U ( i K + 



P + 1 «2 



and u ni m and the time evolution ([22 
is given by 



(23) 



1 



, C f p-i«+ \ / <?-'«+ ) 

T " 2 K+ \^p+iK+ y ((J+i K+ J 



By proper choices of the complex constants C and k + , where k + is defined in the 
upper half plane, and of the real constants p and q one can always render the solution 
U23J) real. 

The two soliton solution is given by 



Pi,: 



E 



P+i 



9+i «j 



J=l 1 



I- 1 k j 
q+i k\ 



with V12 given above. 

4. Continuous symmetries of the lpKdV equation 

Lie symmetries of the lpKdV equation ([2]) are given by those continuous trans- 
formations which leave the equation invariant. From the infinitesimal point of view 
they are obtained by requiring the infinitesimal invariant condition 



(24) 

where 
(25) 



prX„, T 
Mr, 



= 0, 



s 
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By pr X n ^ m we mean the prolongation of the infinitesimal generator X n ^ m to all 
points appearing in D = 0. 

If F n ,m — F n _ rn (u n _ rn ) then we get point symmetries and the procedure to get them 
from eq. (|24p is purely algorithmic [14] . Generalized symmetry are obtained when 
Fn,m = F n , m (u n , m , u n ±i. m , u n ,m±i, ■ • ■)• I n the case of nonlinear discrete equations, 
the Lie point symmetries are not very common, but, if the equation is integrable and 
there exists a Lax pair, it is possible to construct an infinite family of generalized 
symmetries. 

In correspondence with the infinitesimal generator (|25[) we can in principle con- 
struct a group transformation by integrating the initial boundary problem 

(26) du n ,rn(£) _ F n m (u n m (e),u n ±i im (e),u nm ±i(e), . . .), 

as 

with u n _ m (e = 0) = u ntm , where e G K is the continuous Lie group parameter. This 
can be done effectively only in the case of point symmetries as in the generalized 
case we have a differential-difference equation for which we cannot find the solution 
for a generic initial data, but, at most, we can find some particular solutions. Eq. 
(f2~4")l is equivalent to the request that the e-derivative of the equation ID> = 0, written 
for u niTn (e), is identically satisfied when the e-evolution of u niTn (e) is given by eq. 
(|26|) . This is also equivalent to say that the flows (in the group parameter space) 
given by eq. (|26p are compatible or commute with U = 0. 

By means of a symmetry transformation we can construct an e-dependent class 
of solutions of eq. @, given by the functions u„ im (e)- We can associate with this 
class of solutions the corresponding solution tf> n of the spectral problem (|9a[) and eq. 
(|26p can be found among those nonlinear evolution equations associated with the 
spectral problem (|9a|) which are commuting with eq. ([2]). 

4.1. Lie point symmetries. As we do not want to change the basic lattice struc- 
ture of our space, we just consider transformations G e acting on the domain of the 
dependent variables: 

;e), £ e R. 

The infinitesimal generator of the group action of G e on u n _ rn is the vector field 
where 



as 



E=0 

There is a one-to-one correspondence between connected groups of transformations 
and their associated infinitesimal generators [20] . The group action of G e can be 
reconstructed from the flow of the infinitesimal vector field I„ jtI1 either by exponen- 
tiation: 

^n,m (^n,m i — exp(£ X n ra ) 1ln,m 3 

or by solving the differential boundary value problem 

(27) du n} rn(£) = (f, (u nm (e)), 

as 
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with the initial condition u n ^ m (s = 0) = u n ,m- Eq. (|27p is just a subcase of eq. ()26|) 
when F nyTn depends just on u„ )TO . The prolongation of the infinitesimal action of G s 
is given by 

(28) prX n («n 

+ ^ro+lC^n.ro+l) ^Mn, m +i + ^n+l,m+lKi+l,m+l) <9 M „ + lm + 1 . 

Applying the prolongation (|28|) to eq. ((2|) we get 

(29) prX n , m B 

D=0 



— [^n,m(Vro) — <^n+l,m+l(Wn+l,m.+ l)] (p ~ 9 + Un,m+1 — "n+l,?n) + 



+ (p + <7 — U„+l,m+l + Vm) [0n,m+l( u n,m+l) _ <^n+l,m(w n +l,m)j 



= 0. 



j=o 



Due to eq. (|2|) just three of the four different fields appearing in eq. ((29|l are 
independent and, without loss of generality, we choose them to be u n , m , Un,m+i and 

From eq. p9j) . by differentiation with respect to u nm , we obtain, when p =/= q, 



(30) 
so that 
(31) 



fn+l,m+l 



<2m„ . 



du n 



dui 







+ l,m + l 



where m and /C„ m are two integration constants depending just upon n and m. 
Inserting <f> n , m given in eq. (|3 1 [) into eq. (|30j) we find that m = /C„ +1 m+1 , while 
if we insert it into eq. (12^)) we get an explicit equation in w njTO , u n , m +i and u n +i i7n . 
Then the various powers of it n ,-mi M n,m+i and u n+ i >m give us a set of coupled linear 
partial difference equations for JC„ m and m , whose solutions are: 

1C°n,m =71 + (-l)" +m [72 - (?m +pn) 73 ], 

- (-l)" +m 73, 

where 71,72,73 are arbitrary constants. Therefore the lpKdV equation admits a 
three dimensional group G e of Lie point symmetries, whose infinitesimal generators 
are 

(33a) Xl, m = d u n>m , 

(33b) Xlm = (-l) n+m du n>m , 

(33c) X* >m = (-l) n+m [u n , m -(pn + qm)]d Un:m . 

The commutation relations between the generators X^ m5 m , Jf^J m are: 



n.m ' n.m 



0. 



yl y"3 

n,m> n,m 



n.m > n.m 



XL 



so that they span a solvable Lie algebra. 
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The finite transformation generated by the symmetry generators (|33| is given by: 

(34) ^71,711 — 6 ^ ^ ^n,m "i~ 

+ ^'(-1)"+™ 1 [£l+£2 { ~ ir+m ~ ^ (- 1 T +m (P n + 9™)]' 

where £i, £2, £3 are the group parameters associated respectively with the symmetry 
generators X^ m , X^ m) Xf l m . 

When p — q eq. ([2]) simplifies to the product of two linear discrete wave equations: 

(u n ,m — Wn+l,m+l + 2 p) (u n +l, m ~ Un,m+l) = 0, 

and the symmetry group becomes infinite dimensional. 

Finally, let us recall that eq. ([2]) admits also some discrete symmetries: 

(1) The exchange of n and m together with the exchange of p and q leaves the 
lpKdV equation invariant. 

(2) The exchange of n and n + 1 together with the transformation of p into — p 
leaves the lpKdV equation invariant. 

(3) The exchange of m and m + 1 together with the transformation of q into — q 
leaves the lpKdV equation invariant. 

4.2. Generalized symmetries. A generalized symmetry is obtained when the 
function F n m appearing in eq. (|26p depends effectively on u n ^ rn evaluated in some 
lattice points. A way to obtain it is to look at those differential-difference equations 
(|26| associated with the spectral problem (f9"a|) which are compatible with eq. ^ . If 
the nonlinear differential-difference equation (f2"6"| is associated with a spectral prob- 
lem, the proof that it is a symmetry of eq. (|26[) can be greatly simplified. In fact in 
this case the spectral problem establishes a one-to-one correspondence between eq. 
(|26|) and the spectral data and thus one can prove the commutativity of the flows 
in the space of the spectral data, where the equations are linear. 

With this aim, in the following, we construct the hierarchies of nonlinear equations 
associated with the spectral problem (f9"a[) of the lpKdV equation To do so, let 
us rewrite eq. (|9a[) as an operator equation 

(35) L n ijj n = A ip n , L n = E* - {r\ n ^ m + 2p)E n , 

where E^ s ip n = ipn±s, s E N, are the finite shift operators in the n variable the 
eigenvalue A e C is defined in eq. (|12a[) and r\ ntTa = it IljTO — u n+2 , m is an asymptoti- 
cally vanishing potential. Starting from eq. (|35p we can apply the Lax technique, to 
obtain the recursion operator C n and the hierarchy of nonlinear evolution equations 
associated with it, as it has been done in ref. [2]. By the Lax technique we mean 
that constructive procedure introduced by Bruschi and Ragnisco [4 in 1980 which 
consists in finding all those operators M n which are compatible with eq. (|35|) such 
that dtp n /d£ = —M n %j) n . According to the dependence of A on the group parameter 
£ we can distinguish between two different operator equations for the compatible L n 
and M n ; 

(1) if dX/de = 0, then 

(36) — ^- = [L n ,M n ], (isospectral hierarchy); 



CONTINUOUS SYMMETRIES OF THE LPKDV EQUATION 



11 



(2) if dX/de ^ 0, i.e. dX/de = /(A, e), then 

(37) ^p^- = [L n ,M n ] + f(L n ,e), (non — isospectral hierarchy). 

de 

Each one of the two cases will provide symmetries for the lpKdV equation ^ . 

4.2.1. The isospectral hierarchy. The isospectral hierarchy of equations is obtained 
looking at those differential-difference equations whose Lax equation is given by 
eq. (f3T>)) . As we are considering differential-difference equations, M n is a suitable 
operator expressed in terms of shift operators {E^ s }. Taking into account eq. ([35]) , 
the left hand side of eq. (136]) is given by — (dr] nym (e) / de) E n , i.e. it is an operator of 
the form V n (e) E n , where V n (e) is a scalar function. The Lax technique described 
in [4] allows us to write down the following relation: 

(38) V n =C n V n . + V<?\ 

where V n = — (dr) n _ rn (e)/de), V^ -* is a given function of r\ n .m. and of a certain number 
of arbitrary constants. In this case K (0) reads [2] 



(-l) n ai+oa lq n , m + 2p+2j2(- 1 ) k Vk, 



fc=i 



+ a 3 (-l)" I g„, m + 4pn + 2^% >m 

V fc=i 

with a±, a2, 0.3 constants and with q n ^ m defined in eq. (|lip . The recursion operator 
C n is defined by [5] 

(39) C n = -q n , m AH (A1+))- 1 9n , m ^A^))" 1 (AW)" 1 , 

where a1 ±} = ± 1. The inverses of Al ±) can be easily computed in the space of 
functions bounded at infinity and are given by 

00 00 

k=0 k=0 

Choosing V n = 0, from eq. ([38]) a first equation is given by V n — Vn . Iterating this 



procedure we obtain the following set of equations 

(40) ^h^ = c k n v^\ ken. 

dek 

Eqs. ([40]) involve always at least a summation and thus are not local. Following 2J 
we can obtain a set of local equations if we consider in eq. (|40|) powers of the inverse 
of the recurrence operator C n . The explicit form of the inverse operator C~ l can be 
obtained from eq. ([55]) and reads 



= -E-* Al") AW J- (Al-))- 1 AW J_. 

Qn.m Qn : m 

Since A„ = 0, where a is an arbitrary complex constant, we get that (A„ ) -1 = 
a (we can set a = 1 without loss of generality). So applying £~ x to we get a local 
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formula for drj n _ rn /dek- By a recursive application we find 

(41) = c-k-i o. 

ask 

Eqs. (|4ip turn out to be local and they form the so called inverse hierarchy. Taking 
into account the form of C^ 1 we see that with any equation given in formula (|4ip 
we can associate an equation for the field u n>m given by 

(42) ^^c- k ^—+(3 k , 
where 

= -^- 1 -^-A(+)(A(-))- 1 J-A(+)A(-). 

Qn , m y n , m 

Here the ft's are some integration constants to be defined in such a way that u n ^ m , 
asymptotically bounded, is a compatible solution of eq. (|42p . In correspondence 
with eqs. (|42|) we obtain the following evolution of the reflection coefficient 

k 

i n -+- 1 k. \ 

1 - 



(43) '"' 



P 



p — IK 



Taking into account eqs. (pi!?)) and (f4"3")) it follows that 

(44) ^dR^^dR^ yk( _ K 

dEk dek 

where the m-shift operator is defined by E^f f 1h , n = f n ,m±s, sGN, i.e. eq. (HH) is 
a symmetry of the lpKdV equation @. Moreover, we also have 

(45) » U£B , 



dsn V. dek J dek \ deh 
i.e. the symmetries commute among themselves. As eq. ()33a[) is a symmetry of 
the lpKdV equation @ and the symmetries span a Lie algebra, also eq. (l4"2j) with 
/3fc = is a symmetry. However, when ft = the spectral transform is not well 
defined. 

In correspondence with any of the eqs. (|42j) with ft = we can construct an 
infinitesimal symmetry generator 

(46) X^=Z~ k ^— d Unim . 

Qn—l,rn 

Thanks to eqs. (|4"4"H4"5)) the infinitesimal generators (|4T)]) satisfy the symmetry con- 
dition flMJ) and [ X ( n ] , X^ } = V k, h G N. 

For clarity of exposition let us write down the first three equations of the inverse 
hierarchy (for simplicity in terms of q n ,m given by eq. (JXTJ) ) . 
For k = we get 



dUyi m 1 1 



de q n -i,m 2p 

For k = 1 we have 

dUn.m 1/1 1 \ 1 
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For k — 2 we get 



2 

1n-l,m 



+ 



1 



Qn— 2,m \ Qn,m Qn — 1 , m. Qn— l,m Qn—2,m 



l 



i 



l 



Qn — 2 , m Qn— 3,m 



i 



+ 



16 p 5 



QVi.m \ yn+l,m Qn 7 m Qn,m yn— l,m yYi — l,m yVi — 2,m , 

Before going over to the non-isospectral symmetries let us notice that: 

(1) As any of the previous equations can be rewritten just in terms of the variable 
q n , the mappings s„ = (2p)/q n and a„ = s n s n _! enable us to get the 
corresponding equations for the discrete KdV and Volterra equations [7J [TB] . 

(2) Due to the discrete symmetry n «-» m, p <-> g, there exists a symmetric 
isospectral hierarchy associated with the spectral problem (jlOp which give 
rise to a new denumerable class of commuting symmetries. These symmetries 
are given by eq. (|42| with q n . m substituted by 



Pn 



2? - tI n ,m+2 + 



and (E 1 *, A^, £„) substituted by (E±, A, 
symmetry generators of this hierarchy are given by 

1 



(±) 



The infinitesimal 



v( k ) = i — k . 



Pn 



which are also commuting among themselves, i.e. [X m ,X : 



(fe) ^(h)- 



Vk,h G 



N. Denoting by Tk the continuous group parameter of this sequence, the 
first two symmetries of this class are thus given by: 



du r . 



1 



dr Q 
dr\ 



Pn.m—l 
1 



P 2 

I 11,771- 



1 

"2? 
1 



Pn,m—2 



l 

4^3' 



(3) It is easy to prove that the two classes of symmetries we have constructed 
commute among themselves, i.e. [X^\ xffl ] = Vfc, h G N, when D = 0. 

4.2.2. The non-isospectral hierarchy. The non-isospectral hierarchy of equations is 
obtained by looking at those equations which are obtained from the Lax equation 
(|37[) . Let us denote with a (instead of e) the continuous group parameter in this 
case. 

The nonlinear differential-difference equations of the non-isospectral hierarchy are 
given by eq. ([38]) with Vn — q n . Notice that V r i°' ) is not asymptotically bounded, so 
that the resulting difference-differential equations obtained by applying the recursion 
operator C n are not bounded in the limit |n| — > oo. However, also in this case we 
can obtain asymptotically bounded equations by considering the inverse operator 
C" 1 . Taking into account that the solution of the equation A„ ' % = j3 gives 
(A„ ' ) ) _1 /3 = Pn + 7, where j3 and 7 are constants with respect to the variation 
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in n (but can depend on a or m), we can construct a well defined non-isospectral 
hierarchy of asymptotically bounded equations, starting from and choosing 

/(A) = 2A- fc + keN - 



In this case we get 



(47) =£f = CK: 



drj n ^ m k f j 1 



->2 



The first equation of the inverse hierarchy (I47[) is obtained for k = 0. The resulting 
equation is a non-autonomous deformation of eq. ([3]). It reads 
d-g n . m _ In - 1 _ 2n + 3 g„, m 

rfcTO <7ri-l,m 9n+l, m P 2 

which in terms of the field u n ,m is given by: 

dlln tn *2 Tl Un^m. 1 2 71 1 



^ da 2p 2 ' q n -i >m 

The higher order equations of the non-isospectral hierarchy (|4"T)) are all nonlocal. In 
correspondence with eq. ([48]) we have the following evolution of the spectral data: 

(49) ^k= ^ = ^ R 

da aa p{p* + K z ) 

with dn/dcTo = n/p . 

We can now look for (n, m)-dependent symmetries working at the level of the 
spectrum, namely by looking at those evolutions of the reflection coefficients in 
correspondence with equations of the hierarchy (|47|) which commute with the discrete 
evolution of the reflection coefficient given by the lpKdV equation ([2]) and which thus 
satisfy cq. (|4"4"|) with Sk substituted by Ok. If we consider eq. (|4"5)) . as it is the only 
local one in the hierarchy (|47|) it is easy to see that 



dR m dR m +i 
d<jQ dao 

To get a symmetry we could add to eq. (f4"5| any equation of the isospcctral hierarchy 
multiplied by an arbitrary function of m and k, as it has been done in the case of 
the Toda Lattice [7]. To do so, let us add an unknown arbitrary function 7 m (re;p, q) 
to the evolution ([4"9f . 

dRm. ( IK \ 

^ m {K;p,q) R r , 



da \ p(p 2 + K 2 ; y 
and let us require that E m {dR m / d/j$) = dR m +\ / 'dao holds. In such a way we find 

2 i k m 
(q^ + k z ) q 

It is easy to prove that it is not possible to write the function 7 m («;;p, q) given by 
eq. (fST)]) as a finite combination of isospectral evolutions (14"5)) . So we are not able to 
construct in this way non-isospectral symmetries for the lpKdV equation. 
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If we are not interested in looking for asymptotic bounded equations we can 
extract from the non-isospectral hierarchy (|47|) the following local equations 



(51a) 
(51b) 



d,U r , 



pn 



for 



for 



dX 

dlUo 

dX 
dioi 



P 
1 

A ' 



du>i q n -i,m 

where the group parameter is now denoted by to. Eqs. ()51|) are not symmetries for 
the lpKdV equation ([2|). It is possible to show that eq. (I51bjl is a master symmetry 
[Hj of eq. ©. In fact, let us define the following infinitesimal generators 



(52) 
(53) 



y(i) 

n 



Qn —l,m 



pn,m—l 



where Y m X ^ is obtained from Yn by applying the discrete symmetry n <-> m, p <-> q. 



As Yri 1 ^ (or equivalently Y m L> ) is not asymptotically bounded, one cannot associate 
with the corresponding equation (|26[) the evolution of a reflection coefficient and 
thus we cannot compute the commutation relations [Xn ,Yn ] for an arbitrary k. 
Let us write down the simplest commutation relations: 

3 



Hi) 



y-(O) y(l 

n 1 -* n 



A similar result could be obtained considering the commutation of X m and Ym . 
As one can see from the above commutation relations the commutations of Xn with 
Y„ and of with Ym provide an alternative constructive tool with respect to 
the use of the recursive operator to find the infinite families of symmetries X n k ^ , xffi 
for the lpKdV equation @. 

4.2.3. Construction of non- autonomous generalized symmetries. We can now con- 
struct some new (n, m)-dependent non-autonomous symmetries of eq. Such a 
construction is provided by the following theorem whose proof is straightforward. 

Theorem 1. Let H(u n ,mi ^n±i,m) Ura,m±i> •••;£>) q) — be an integrable partial dif- 
ference equation invariant under the discrete symmetry n <-> m, p «-> q. Let Z n be 
the differential operator 



such that 
(54) 



pr Z„ 



i' = 



— Qn^m \U>n,mi ^n±l,m , Wn,m±l , •••) Pi q) t 



where g n ,m( u n,mi u n±i,m> w n ,m±i> •■■jPi q) * s invariant under the discrete symmetry 
n <-> m, p <-> q and a is an arbitrary constant. Then 



1 



■ pr Z n 



prZ n 



0. 



5=0 
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where the operator Z m = Z m {u n>m , ti n>m ±i, u n ±i j7n , q,p) d Un m is obtained from 
Z n under n «-> m, p «-> q, so that 



— P 9n,m (^n.m 3 ^n±l,m 3 ^rl.mil 3 ■•■iPi q)i 

=0 



pr Z m D 
being (3 a constant. So 

(55) Z nm = Z n — Z m 

a p 

is a symmetry o/B = 0. 

The above theorem provides a constructive tool to obtain generalized symmetries 
of the form given in eq. ([54"]) for the lpKdV equation © . Such a situation occurs with 
the operators Y„\Ym in eqs. (p32"H53"]) as in this case ag n . m = —1 and g n , m = 1- 
The symmetry + Ym has been considered by Tongas in |23j and the invariance 
condition one gets from it has been derived by Nijhoff and Papageorgiou [TS] from 
a monodromy problem associated with eq. (|8a[) . 

A more complicated symmetry for eq. ^ can be obtained from Theorem 1 by 
defining an operator Z^ , depending parametrically on w G M, with 

7 (™) ■ np w p w -q w ( 1 

56 Z\ > = — -2 2T pn-- u n>m 

Sn-l.m 2 (p^ - q*) \ 2 

In such a case we find: 

a£ m(p W - <f) + 4pqx n , m {p w - 1 + q 1 "- 1 ) + Ap 2 q 2 (p w - 2 - q w ~ 2 ) 



4 (p + q) (x 

II . II! -p + q) 



where x n>m = u n+ i, m - u n>m+ i. The operator Z n = Z^d Un m , with Z^ given 
by eq. ([56]) . satisfies the condition ([54]) and thus the operator f|55[) defines a one- 
parameter symmetry for the lpKdV equation For w = in eq. (|56p we recover 
the simmetry generated by In 1 ** + Ym^ ■ The simmetry obtained by setting w = 1 in 



eq. ([56]) can be also found in [23] . 

Let us give here two other examples of functions Zn satisfying condition ([54]) : 

(57) z {w) = n (P W + 1 W ) _ 

q n -l,m p w -q wl 

(58) = n (P q } W _ 

q n -l,m P w -q w ' 

with Ul6l\ {0}. 

Using Theorem 1 we can also show that it is not possible to construct a symmetry 
using eq. (|51a.|) . In fact, defining the differential operators 

Y n 0) = ( U n,™ ~ P Tl) Un m , Y^ = (u n<m - q m) d Un m , 

we immediately get that condition (|54|) is not satisfied as 

_ qxl, m + 2 P(p-q) x n , m - 2p 2 (p-q) 

Q! Qn.ra j 
>£n,m P ~r Q 

is not invariant under the exchange n <-> m and p ^ q. 
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4.2.4. Symmetry reduction on the lattice. Let D = be a discrete equation and 

•)du n ,m one °f its symmetry generators. A 
solution u ntm of D = is a solution invariant under X n _ m if it satisfies the equation 

We now use the notion of invariant solutions of the lpKdV equation ([2|) to con- 
struct a particular solution. Our results are a slight generalization of the ones 
contained in [19l [23] . 

Hereafter, carrying out the simplifying transformation 



(59) \/p~+~q +pn + qm, 
we shall write the lpKdV equation ([2|) as 

(60) (u n+ i <m - U n , m +i) (w„+l,m+l - Un,m) = P ~ 1- 
Let us consider the symmetry generator X n _ rn = y/p + q Z^m + cX\ m , 

where ZnX = Zn + Zm\ with given in eq. ([56"]) and X\ m given in eq. (|33a|) . 

Under the map ([55]) , the differential operator X n m reads 

(61) 

9 _( np w mq w 1 p w - q w \ 

V^n+l.m ^n— l,m ^n,m+l ^n,m— 1 ^ P ^ / 

Solutions invariant under the symmetry generator (|6ip are solutions of eq. (|60p 
subject to the constraint 

(62) + 



^n+l.m ^n,m Vri.m Vn 



^n+l,m ^n — l,m V>n,m+1 ^n,m— 1 2 p ^ 

Following [13] we define a„ !m = m„+i, to - u n -i,m.) &n,m = Un,m+i - «n,m-i and 
2M,m = u n+ i, m+ i - u n , m . Using eq. ([10]) we get 

(5 

^n.m — l,m "T j 

with 5 = p — q. On the other hand, we have 
(63) 

Eliminating the variable b n ^ m between eqs. (|62H63]) we obtain the following second- 
order difference system: 

p w (n+l)5 p w n8 w . .,. w 5 

— '— + — - r = p w (n + 1) + q w rn - c— + 

VnVn+1+0 VnVn-l+O ' y n 

(64a) + cy n H — — [y n u n + 6 

(64b) u n+ i - u n _! = j/ n _i H , 

where we have dropped the index m in the dependent variables as it enters just 
parametrically. Eqs. (|64aH64bf are a system extension of the alternate discrete 
Painleve II equation [5] Q]|] . 



18 



DECIO LEVI AND MATTEO PETRERA 



A similar procedure can be performed considering the operators (|57l l58f . In these 
cases one obtains, as a reduction of eq. (|60j) . an equation that is equivalent to the 
alternate discrete Painleve II equation for any value of w. 

5. Concluding remarks 

In this paper we have studied in detail the lpKdV equation and its symmetries. 
To do so in a complete way we started from its associated spectral problem. The 
study of the direct and inverse problem allow us to provide the soliton solutions of 
the nonlinear lattice equation. As has been shown in [TTl |T5] the associated spectral 
problem of the lpKdV equation is given by the asymmetric discrete Schrodinger 
spectral problem considered by Boiti et al. [3]. A class of generalized symmetries of 
the lpKdV equation is obtained by considering its associated isospectral nonlinear 
evolution equations. 

It is worthwhile to notice here that the same transformation that exists between 
the Volterra equation ([7]) and the discrete KdV equation $Q can be written between 
the corresponding spectral problems. So part of the results we presented could be 
obtained from the discrete Schrodinger spectral problem associated with the Volterra 
equation (J7J). The eq. |6]) is associated with the spectral problem (f5]) which, in the 
variables of the discrete KdV equation, reads 

2p 

(65) tp n+ 2 = tpn+1 + A^„, 

while the spectral problem for eq. ([7]) is 

(66) 4> n -i + a n 4> n+ i = Lh(/> n . 

By up-shifting eq. (|66p by one, defining a n = s n s n _i and considering the trasfor- 
mation <j> n = f n ([i) <? n ( s n> s n±i> •••) we are & ble to rewrite eq. (|66p as eq. (p5|) . 
To do so we choose /„ = fi n with A = — 4p 2 //i 2 and g n given by the solution of 
the recursion relation g n +i — <?n/(2ps„). From this mapping, but in particular 
from the relation between A and [i, we get that the positive Volterra hierarchy [7j 
corresponds to a negative hierarchy in the asymmetric discrete Schrodinger spectral 
problem (|65p . In this way the Miura transformation between the potential a n , m of 
the Volterra hierarchy and field u n ,m of the lpKdV equation is easily obtained and 
is given by 

4p2 

n ' m (2p- U n+2 ,m + u n ^ m ){2p- U n+ \ >m + U„_l >Tn ) 
From the two equivalent components of the Lax pair we get two infinite hierarchies 
of integrable symmetries by looking at its isospectral deformations. From the non- 
isospectral deformations we get an infinite set of master symmetries. Apart from the 
two classes of integrable symmetries, associated with the two equivalent members 
of the Lax pair, we have obtained an infinite set of symmetries starting from the 
master symmetries. It is an open problem to prove whether the symmetries Z^ w > are 
integrable, possibly by providing their Lax pair and to compute their commutation 
relations. In any case we can use them to provide solutions by symmetry reducing 
the lpKdV. 

In conclusion we have shown that the lpKdV equation has a lot of symmetries. 
It is left to future work the proof if all are independent or if, the two sequences 
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associated with the spectral problems in n and m respectively are dependent under 
the equation. 

It is also to be clarified the role of the Z^ symmetries, related to the master 
symmetry. It is an open problem to see whether the other discrete symmetries might 
give rise to new continuous symmetries. We leave to a future work the use of the 
symmetries here obtained to check the integrability of the (local) dNLS equations 
derived in [§1 QUI Q3] . 

After we finished this article a referee for our recent paper [13] pointed out to us 
the article by Rasin and Hydon [21] submitted to Stud. Appl. Math, and available 
in the personal web page of the authors. In this article the authors construct the 
symmetries for all quad-graph equations contained in the Adler-Bobenko-Suris list 
P]. Let us notice here that in the present paper we have just considered the Hi 
equation but with much more details. In particular we already answer in this case 
to some of the questions in their conclusions. Moreover we have obtained five point 
symmetries that they have not been able to find. 
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